Dark energy and cosmic microwave background bispectrum 
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We compute the cosmic microwave background bispectrum arising from the cross-correlation of 
primordial, lensing and Rees-Sciama signals. The amplitude of the bispectrum signal is sensitive 
to gthe matter density parameter, flo, and the equation of state of the dark energy, which we 
parameterize by wq. We conclude that the dataset of the Atacama Cosmology Telescope, combined 
with MAP 2-year data or the Planck data set alone will allow us to break the degeneracy between 
Qo and wq that arises from the analysis of CMB power spectrum. In particular a joint measurement 
of Qo and wq with 10% and 30% error on the two parameters respectively, at the 90% confidence 
level can realistically be achieved. 
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I. INTRODUCTION 



Recent results suggest that the Universe is flat and 
dominated by a negative pressure (dark energy) compo- 
nent (e.g., [1] and references therein) which can be char- 
acterized by its equation of state parameter p/p — wq. 
A constant vacuum energy (i.e. a cosmological con- 
stant) has WQ = —1, while a generic wq < dark 
energy component is refereed to as "quintessence" (e.g., 
[2]). Analysis of the power spectrum alone of forthcom- 
ing Cosmic Microwave Background (CMB) experiments 
will still present some degeneracy between wq and fio 
(after marginalization over other cosmological parame- 
ters) (e.g., [3]). 

Here we consider the constraints on wq-^Iq that can be 
obtained from the analysis of the small scale CMB bis- 
pectrum for two experimental settings: the combination 
of MAP [4] and Atacama Cosmology Telescope (ACT; 
[5]), and the Planck surveyor satellite [6,7]. 

Under the assumption of Gaussian initial conditions, 
the cross-correlation of primordial, gravitational lensing 
and Rees-Sciama (RS) [8] signals is the dominant contri- 
bution to the CMB bispectrum after that of the Sunyaev- 
Zeldovich (SZ) effect [9] and of point sources. Since each 
term in the bispectrum has a different function form, we 
can separate these signals without major loss of informa- 
tion [e.g., 10]. 

This paper is organized as follows. In section 2 we 
present the expression of the primordial-lensing-RS cross- 
correlation bispectrum. In section 3 we forecast the error 
on the joint determination of fio and wq from the bispec- 
trum analysis. In section 4 we conclude that it is possible 
to break the degeneracy between flo and wq that arises 
from the CMB power spectrum analysis alone. 
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II. THE 

PRIMORDIAL-LENSING-REES-SCIAMA CROSS 
CORRELATION BISPECTRUM 

We wish to compute the effect on the CMB bispec- 
trum of the coupling between the RS and the weak lens- 
ing. The weak lensing effects allow us to probe the low 
redshift Universe through the lensing imprint that fore- 
ground structures leave on the primordial CMB signal. 

Neglecting galactic contamination, the CMB tempera- 
ture at any position in the sky 7 can be expanded as: 
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where the first term is the primordial signal, the second 
term is the gravitational lensing effect, the third term is 
the RS contribution and the last term denotes the SZ 
effect. The frequency dependence of the SZ term makes 
it possible to separate out this contribution and we will 
therefore ignore this term in what follows. Point sources 
contribution to the bispectrum signal can be separated 
out without major loss of information [10]. Other con- 
tributions to AT/T such as e.g., the Ostriker-Vishniac 
effect [11] will have zero or sub-dominant contribution to 
the CMB bispectrum. 

The primordial contribution can be expressed as: 
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where g denotes the radiation transfer functional, $ de- 
notes the Fourier transform of the gravitational potential 
perturbation $, and denotes the conformal distance to 
the last scattering surface. The lensing potential 0(7) is 
the projection of the gravitational potential: 
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The RS effect arises from a combination of two ef- 
fects: the late-time decay of gravitational potential fluc- 
tuations in a non-Einstein-de Sitter Universe — strictly 



2 



called Integrated Sachs- Wolfe effect [12] — and the non- 
linear growth of density fluctuations [8] along the photon 
path. The third term in Eq. (|l|) can be expressed as 
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The CMB bispectrum is defined as: 
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where a™ are the coefficients of the spherical harmonics 
expansion of the observed CMB temperature fluctuation: 
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) is the Wigner three-J symbol, and the last 
equality results from symmetry reasons (e.g., [13,14]). 
By applying Eq.(|) to Eq.(|lj) we obtain (cf. [15,10]), 



and Qi are the normalized densities of the various energy 
components of the Universe. The exponents depend on 
how each component density varies with the expansion 
of the Universe, pi oc a"', where a is the scale factor 
and Hi — i{l + Wi) , and consequently Wi is the equation 
of state parameter for the i component. For example for 
matter n — 3, for a cosmological constant w — —l,n — 0. 



A. Computation of Q{£) 

We compute Q{£) numerically for different combina- 
tions of Oq and wq for COBE-normalized models. We 
assume a flat Universe and set fth^ = 0.17. This is justi- 
fied by the fact that, from MAP 2-year data, this quantity 
should be known to better than 5% accuracy [17]. 

The gravitational potential power spectrum at any 
given redshift z is related to the density power spectrum 
(P) via: 
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where TL denotes Gaunt integral 
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Following the steps outlined in the Appendix we obtain 



an expression for 2(^3) = (9 
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(cf. [15,16]) 
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where z* denotes the redshift of the last scattering surface 
and P$(fc, z) denotes the power spectrum of the gravita- 
tional potential at redshift z; it has to be evaluated at 
k = l^/r and then derived with respect to z. 

Since we assume a flat Universe, r{z), the conformal 
distance from the observer at 2: = is 
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where 
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In the linear regime 
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and for the bispectrum we obtain 

_ ^^1^2^3 ^l(^l + 1) - ^2(^2 + 1) + ^3(^3 + 1) 

- '^{■11213 2 

X Ci;" (9™^' a^^^™3) 5 permutations . (9) 



P^'^iKz) = Ak-^Tikf (iW il±ffl) (14) 

\g(z*) (1 + 2) / 

where T{k) denotes the matter transfer function, g{z) 
the correction to the growth factor due to the presence of 
dark energy, A is the amplitude of the primordial power 
spectrum (see [18]), and Ug denotes the primordial spec- 
tral slope. 

In the case of wq = —1, (i.e. for a cosmological con- 
stant) the expression for g{z) is: 
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otherwise there is a correction factor [19] (— wq)*, where 

t = -(0.255-|- 0.305u;q-K 0.0027/wq)(1-17;,) 
- (0.366 + 0.266WQ- 0.07/^(2) In r?^. (16) 

Here Vl^ is = Qq/ [^q + (1 - f^o)a~^'"''] • We approx- 
imate the transfer function with that of Sugiyama [20]. 
Any corrections for w 7^ — 1 affect only very large scales, 
that do not contribute to the signal we are modeling. 

Since the signal for Q is mostly coming from non-linear 
scales, Eq. (^ is not a good approximation of the actual 
power spectrum. The nonlinear mass power spectrum 
can be obtained for the linear one with the mapping of 
Peacock and Dodds [21] for wq ~ —1 while to generalize 
the mapping to wq 7^ — 1 we use the expression of Ma et 
al. [18]: 
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FIG. 1: Absolute value of Q{1) for two different cosmologies: 
^0 = 0.3, w = -i.O (thick line) and fio = 0.2, w = -0.2 
(thin line). The solid line indicates Q > while dashed line 
indicates Q < 0. If linear theory was a good approximation 
for the evolution of the power spectrum, Q > 0: non-linear 
effects can be important at £ ~ 200. 



where fci = fc/(l + ^2 ^ 2/T:^k^P{k), 



G(x) = (l + ln(l + a;/2)) 



1 + 0.022;'' + cixV^^ 
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, (18) 



ci = 1.08 X 10""*, C2 = 2.1 X 10"^ and qa = 

|^|1.3h|-0.76^(0). 

The sign of i9P$ /dz in the integrand of Eq. (|l^) is de- 
termined by the balance of two competing contributions: 
the decaying of the gravitational potential fluctuations as 
z — > and the amplification due to non-linear growth. 
Both of these are sensitive to the cosmological parame- 
ters, we thus expect Q{£) to be sensitive to wq and Uq. 
In Fig. (|l|) we show the effect of non-linear evolution on 
Q{£). The solid line indicates Q > while dashed line 
indicates Q < 0. If linear theory was a good approxi- 
mation for the evolution of the power spectrum, Q > 0: 
non-linear effects can be important dX (. ^ 200. 



III. A PRIORI ERROR ESTIMATION 

We now wish to evaluate how well forthcoming CMB 
experiments will be able to measure Jig and wq, if we con- 
sider the information enclosed in the primordial-lensing- 
RS bispectrum in addition to the power spectrum. We 
will thus estimate the on a grid of cosmological models 
(cf. [10,15]): 
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where for symmetry we consider ti 
all ts are different, = 2 if two £'s are repeated and 
A'' = 6 if all Ps are equal. In deriving Eq.(19) we have 



used the fact that the sum over m of the square of the 
Wigner three-J symbols is unity. 



Following the statistic introduced by [15] the confi- 
dence region for fio and wq jointly at the 68.3%, 90% 



and 99%, is then assumed to be given by 5x 



=2.3, 4.61 and 9.21 respectively. Of course, the 



relation between Ax^ and confidence levels is strictly cor- 
rect only if our "data" Bg^g^g^ are normally distributed. 
The central limit theorem ensures that for a large number 
of independent data the distribution should asymptoti- 
cally tend towards a Gaussian, nevertheless this assump- 
tion will need to be tested a 'posteriori or a maximum 
likelihood technique will need to be used. 



In computing we make several approximations: first 
of all the expression for Q (Eq.|l^ uses the small angle 
approximation. For the purposes of error estimation this 
approximation is good for i > 10. However most of the 
signal comes from the coupling of very large scales (small 
i) to very small scales (large £). Since the exact ex- 
pression is computationally expensive, we consider only 
i > 12 in our x^ calculation on the grid of cosmological 
models. For the standard ACDM model, using the exact 
expression for Q at small £'s, we obtain that the x^ is 
amplified by a factor 2if2<£< 12 are also included. 
We can thus infer that this approximately applies to all 
flo, Wq combinations ad that the cosmological parame- 
ters determination can be improved consequently, if we 
include in our analysis all £ > 2. 

In Eqs. ( p^ ) and ( p^ ) we approximate the power spec- 
trum (Cg) as composed by three contributions: primor- 
dial (Cf ), Ostriker-Vishniac (Cf ^) and noise (Cf ). C[ 
is obtained using CMBFAST [22] up to ^ = 1500 (with 
parameters ilo — 0.3, A = 0.7, = 1, and, conserva- 
tively [ p6| , T = 0) and is approximated by a power law 
for e > 1500. The OV contribution (AT°^) is conserva- 
tively taken to be 2.8 x 10~^ y/2TT/P; i.e. ^ becomes 
important only at ^ > 3000. For the noise calculation 
we assume that the experimental beam is gaussian with 
width CTb ~ Ofwhm f'^'-i, where OpwHM is the beam full 
width at half maximum. Following Knox [23] we have 
that cf — exp(^^cr^)5, where 5* is the instrument sen- 
sitivity i.e. the noise variance per pixel times the pixel 
solid angle in steradiants. For the noise contribution from 
many independent channels (as for Planck case for exam- 
ple) we use (Cf )-i = T.A^^Mi'^))'^ ■ We consider two 
different experimental settings. One with Planck speci- 
fications the other with ACT for £ > 200 and MAP for 
£ < 200 specifications. ACT will map about 1/100 of 
the sky with 2fiK per pixel noise and experimental beam 
with OpwHM — 1-7'. Details about Planck specifications 
can be found e.g., [7]. In practice, for the combination of 
ACT and MAP datasets, useful signal can be extracted 
up to ^ ~ 9000 while for Planck up to ^ ~ 2000. 
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FIG. 2: Confidence levels in the ^o-wq plane. The dashed 
and dotted contours show the degeneracy in the plane for the 
2-year MAP power spectrum data (see text for details). The 
other lines show the expected confidence levels from the bis- 
pectrum analysis described in the text applied to MAP and 
ACT data for different fiducial models (indicated by the di- 
amond). In particular solid lines show the 68.3% (labelled 
by B), 90% (C) and 99% (D) confidence region for Q.q and 
wq jointly considering only ^ > 13 and the small angle ap- 
proximation. We estimated that by considering £ > 2 the 
constraints on cosmological parameters are improved as fol- 
lows: the 68.3% confidence level region is indicated by the 
dot-dashed line labeled by A, the 90% confidence level corre- 
spond to the line labeled by B and the 99% correspond to the 
line labelled by C. Note that for low-f2o fiducial model the 
constraint are much more stringent than for high-f2o model. 



A. Breaking the degeneracy 
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FIG. 3: 90% confidence level joint constraints on fio, wq from 
MAP 2-year data and ACT for two fiducial models (indicated 
by the diamond) . Light gray shaded area is excluded by MAP 
power spectrum analysis alone, while dark gray shaded area 
is excluded by MAP-I-ACT small scale bispectrum consider- 
ing only I > 12. We extrapolate that the area filled with 
pattern can be excluded by considering also 2 < ^ < 12 (see 
text for details). Similar constraints can be obtained from an 
experiment with the specifications of the Planck mission. 



(D) confidence region for fip and wq jointly. To obtain 
these contours the has been computed conservatively 
considering only £ > 12 and the small angle approxima- 
tion. However the constraints on cosmological parame- 
ters can be improved by considering all £ > 2. In this 
case the 68.3% confidence level region is indicated by the 
dot-dashed line labeled by A, the 90% confidence level 
correspond to the line labeled by B and the 99% corre- 
spond to the line labelled by C. Note that for low-fio 
fiducial model the constraint are much more stringent 
than for high-f2o model. 

Similar constraints can be obtained from an experi- 
ment with the specifications of the Planck satellite. 



While observations of the microwave background fluc- 
tuations are sensitive probes of cosmological parameters, 
there are significant parameter degeneracies. In a flat 
universe, the position of the first acoustic peak depends 
primarily on the angular diameter distance to the surface 
of last scatter. In a universe with dark energy, this dis- 
tance is a function of f^o and wq. We have explored the 
degeneracy by simulating a Monte Carlo Markov chain 
analysis of the MAP 2 year data. In the analysis, we have 
assumed that the MAP data is limited by the statistical 
errors and applied the Monte Carlo Markov Method de- 
veloped in [24]. In our analysis, we have fit the data with 
a seven parameter model (power spectrum amplitude, 
power spectrum slope, baryon density, matter density, 
Hubble constant, rcionization rcdshift and wq). 

In Fig. 1^ we show confidence levels in the VIq-wq plane 
for different fiducial models. The two dashed contours 
and the dotted contours show respectively the 99%, 90% 
and 68.3% from the 2-year MAP power spectrum data. 
The other lines show the expected confidence levels from 
the bispectrum analysis for MAP and ACT data. Solid 
lines show the 68.3% (labelled by B), 90% (C) and 99% 



IV. CONCLUSIONS 

We have computed the effect on the CMB bispectrum 
of the coupling between Rees-Sciama, gravitational lens- 
ing, and primordial signal. This signal is determined by 
the balance of two competing contributions along the line 
of sight: the decaying gravitational potential fluctuations 
and the amplification due to linear gravity. Both of these 
effects, and thus the bispectrum itself, depend on wq and 
r^Q. Since most of the bispectrum signal comes from the 
coupling of large scales (low €} with small scales (large €} 
we have examined two experimental settings that allow 
to accurately measure CMB temperature fluctuations at 
low and high ^'s: one is a combination of MAP 2-year 
data with ACT CMB maps and the other has the spec- 
ifications of the Planck surveyor. As shown in Fig. (||) 
we conclude that we can realistically achieve an error of 
about 10% on fig and 30% on wq at the 90% joint con- 
fidence level, by combining the constraints from CMB 
power spectrum and bispectrum. 

It is however important to bear some caveats in mind. 
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In general wq might be time dependent. The CMB power 
spectrum will give constraints on some " weighted mean" 
of wq{z). The analysis presented here constraints a dif- 
ferent weighted mean of 'wq{z), where most of the weight 
comes from z ~ 1. This method has to be interpreted as 
a first order approximation to detect wq 

We have also assumed that the CMB primordial signal 
is gaussian and that other foregrounds contributions to 
the bispectrum (e.g., dust, point sources, SZ effect) can 
be subtracted out. While the SZ and point sources con- 
tributions can be accurately subtracted out (e.g., [10]), 
dust should be negligible above the galactic plane and 
accurate dust templates are available [25], the presence 
of a primordial non-gaussian signal might invalidate the 
results. 
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V. APPENDIX 

The derivation of Eq. ( |lo|) is conceptually similar to 
that of Spergel & Goldberg [15] for the Integrated Sachs- 
Wolfe effect, but is complicated by the fact that the non- 
linear evolution of the power spectrum cannot be factor- 
ized in a time dependent and a scale dependent functions. 
We start from: 



X y'dr<i>^^(r,72r)y;™^(7i)r,7(72)) (20) 

where the dot denotes d/dr. Writing $ in terms of its 
Fourier transform $ and expanding the exponential as 
exp(fc.7r) = {AT:)Y.i'rn'^''Yr' {lk)Yir' mt{kr), we 
obtain: 

X P$(fc; r, ry+'"n. {kr)Y;r' {%)Ylj' (71) 

xr;,™" {%)YiJ}" (72)r™^ {^,)y^^ (72) (21) 

where P$(fc;T, r) is defined through 

(i>(fc,r)$(fc',r)) = 2(27r)3p<i,(fc;T,r)5^(k + k') (22) 

and 5^ denotes the Dirac delta function. In principle 
(21) has an extra term which vanishes at high £. Us- 
ing the fact that d^k = k^dkd% and the orthogonality 
relations of spherical harmonics, we obtain J d^id^2d'^k 
yields -Sp^^d^^^,Sp,^J^^^„Sp,g,6^„_^, , where 6^^ de- 
notes the Kronecker delta. Finally using the approxima- 
tion: 

j dkk'^f{k)jt{kr)ki>,{kT)~ f{i'/r)TT/2r^5{r-T) (23) 

and performing the remaining integral in dr, we obtain 
Eq. 0. 
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